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Abstract 

We map the phase diagram of gauge theories of fundamental interactions in the flavor- 
temperature plane using chiral perturbation theory to estimate the relation between the pion 
decaying constant and the critical temperature above which chiral symmetry is restored. We then 
investigate the impact of our results on models of dynamical electroweak symmetry breaking and 
therefore on the electroweak early universe phase transition. 
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I. EXTREME TECHNICOLOR 

Arguably technicolor extensions of the standard model (HE] constitute a natural mech- 
anism able to explain and drive the breaking of the electroweak symmetry (see the recent 
review [3]). 

As for ordinary Quantum Chromodynamics we can consider technicolor models in 
extreme conditions of temperature and matter density ||4]-[9]]. These regimes are particu- 
larly interesting for the early universe dynamics and phenomenology |6r|8]|. For example 
the technicolor temperature driven chiral phase transition is directly mapped into the 
electroweak phase transition. Applications to the detection of gravitational waves is very 
interesting and can potentially distinguish different scenarios of dynamical electroweak 
symmetry breaking |10f. 

Adding a nonzero matter density is also very interesting for investigating asymmetric 
type dark matter emerging naturally as a technicolor interacting massive particle. The 
early models dealt with scaled up version of the ordinary baryons, i.e. technibaryon IfTTT — 
[T3| . Recent models of technicolor [3 J led to the introduction of new types of dark matter 
candidates [14-24]. Introducing nonzero chemical potentials allows us to investigate 
the dark matter composition and distribution in the universe. An intriguing point is 
that for the recent relevant models of technicolor [3] the action remains positive upon 
introduction of the chemical potential and therefore the nonperturbative dynamics can 
be studied using first principle lattice computations. This is so since in these models the 
sign problem is absent. An initial study of the nonzero matter density and temperature 
phase diagram, highly relevant for the new dark matter candidates, has been presented 
in |25|. Here it was also shown that there is an intriguing relation between the matter 
density driven chiral phase transition and the deconfining one yielding a very rich phase 
diagram as function of the temperature and matter density. 

In this work we start mapping the phase diagram of gauge theories of fundamental 
interactions in the flavor-temperature plane using chiral perturbation theory to estimate 
the relation between the pion decaying constant and the critical temperature above which 
chiral symmetry is restored. We then investigate the impact of our results on models 
of dynamical electroweak symmetry breaking and therefore on the electroweak early 
universe phase transition. 
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II. CHIRAL LAGRANGIAN SET UP 

If the number of flavors is sufficiently low we expect chiral symmetry to break for any 
gauge theory of fundamental interactions we consider here. It is, hence, always possible to 
construct a chiral Lagrangian in terms of the Goldstone bosons of the theory and compute 
the temperature dependence of the chiral condensate ||26l - (32f . At the leading order in the 
chiral expansion the only scale in the problem is the pion decay constant. 

A. Single Matter Field Representation 

First, we will discuss the theories with fermions transforming under a single represen- 
tation of the technicolor gauge group. For a theory with 2Nf Weyl fermions, the chiral 
symmetry G is either SU(Nf) L x SU(Nf) R or extended to SU(2Nf), if the representation is 
complex or (pseudo)real, respectively. We expect that the chiral symmetry breaks sponta- 
neously to its maximal diagonal subgroup H. This leads to the following chiral symmetry 
breaking patterns G — > H: 

• SU(Nf) L x SU(Nf) R — > SU(Nf) v for complex representations, 

• SU(2Nf) — > SO(2Nf) for real representations, 

• SU(2Nf) — > Sp(2Nf) for pseudoreal representations. 

Examples of (minimal) technicolor theories ||3l [T4ll23ll33l - l36l with some of the above chiral 
symmetry breaking structures are: 

• Minimal Walking Technicolor (MWT), which has N f = 2 in the adjoint of SU(2) TC 
with expected pattern SU(A) SO(4). 

• Next-to-Minimal Walking Technicolor (NMWT), which has Nf = 2 in the symmetric 
representation of SU(3) TC with the pattern SLT(2) L x SLT(2) R SU(2). 

We shall construct the nonlinear chiral Lagrangians following Refs. [[141 [151 (other kinds 
of effective theories for minimal walking techicolor models have been considered in ||3l[T4l 



IT51I37H4011). In a basis formed by the Weyl fermions \Q\, . . . , Q L N , io 2 (Q*J ,...,io 2 h 
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the chiral condensate is proportional to a 2Nf X 2NV matrix which we denote by: 

a) 



' 1^ 



±1 

with the plus sign for complex or real representations, and minus sign for pseudoreal 
representations. Here the subblocks are Nt x Nf matrices. We denote the unbroken 
generators by S a . They leave E invariant: (S a ) T E + ES" = 0. The broken generators are 
denoted by X'. Explicit expressions for the generators can be found in |tT5ll37||. Let us 
define the element *V of the coset space 



1 / = exp^n i X i 



E, (2) 



where F n is the Goldstone boson decay constant. It transforms non-linearly 

VW^gVWt&g) (3) 

where g is an element of G and h is an element of H. We can embed the electroweak gauge 
group in SLT(4) as done in |T4"ll3~7ll . It is appropriate to introduce the Hermitian, algebra 
valued, Maurer-Cartan one-form 

to, = W f d,V (4) 
From the above transformation properties of *V it is clear that to, transforms as 

co, ^ g)co,h\Z, g) + g)d,h\Z, g) . (5) 

With co, taking values in the algebra of G we can decompose it into a part or, parallel to 
H and a part co^ orthogonal to H 

J, = 2S"Tr [s a co^ , coj; = 2X'Tr [Xco,] . (6) 

Then co, is an element of the algebra of H while co^ of G/H. We find the following 
transformation properties: 

J, -> h(l, g)co\}x\l, g) + h(l, g)d,h f a, g), cof^ h(l, g)^h\l, g) . (7) 

To probe the chiral dynamics of a theory we introduces a small, with respect to the 
dynamically generated scale associated to the chiral condensate, democratic mass matrix 
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oc E for the fermions at the underlying theory level. Since the mass term introduces an 
explicit breaking of the chiral symmetry the pions become massive. The democratic choice 
of the mass leads to equal masses to the Goldstone bosons. At the effective Lagrangian 
level we have: 



mlF 2 



L = F 2 Jr [co^ x ] - -^Tr [E + T r ET + h.c] . 



(8) 



B. Multiple Representations: The Ultra Minimal Technicolor example 



Lagrangians for technicolor models featuring matter transforming under several dis- 
tinct representations of the underlying gauge group may be constructed in a similar 
fashion. The minimal walking model known as Ultra Minimal Technicolor (UMT) is one 
of the phenomenologically relevant examples. It contains two Dirac fermions transform- 
ing according to the fundamental representation and two Weyl fermions belonging to the 
adjoint representations of the technicolor gauge group SU(2) TC . Its nonlinear Lagrangian 
was constructed in [14]]. UMT has one anomaly free U(l) symmetry. We expect the non- 
perturbative dynamics to spontaneously break the global SLT(4) x SU(2) x IT(1) symmetry 
to Sp(4) x SO(2) x Z2 via the formation of the two distinct condensates. Let us order the 
broken generators such that {Xi, ...,X 5 }, {X 6 ,X 7 }, and X 8 are part of the generators of 
SLT(4), SU(2), and U(l), respectively (see [14j for an explicit realization). We define 



(9) 









1 




'or 


E 4 = 


-1 

v. ) 


; e 2 = 


1 



where 



The element *V can be parameterized by 

^) = exp(^X ! )E, 

where 



(10) 



(11) 



i=1 



ITX 



n 8 x 8 



(12) 
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and F n , F n and F n are the related Goldstone boson decay constants. We introduce indepen- 
dent mass terms for the fundamental and adjoint techniquarks. We follow the procedure 
outlined above to construct nonlinear Lagrangians and by noting that the generator X 8 
corresponding to the broken U(l) is not traceless we deduce: 

X = Tr [aa^a^] + bTr [<»i] Tr - Tr [ cE if V T E<V + h.c] , (13) 

The coefficients a = diag ^F^,F„,F^,F^,F^,F^ and b = y - -y - ^ are chosen such that 
the kinetic term is canonically normalized: 

1 8 

X=-^^,nw + ... . (u) 

The remaining coefficient is c = diag (m^F^ / 4, m^F^/ 4, m^F^/ 4, m^F^f 4, fh^F^/ 4, m 2 n F 2 n l4}, 
where m n (m n ) is the mass of the pions emerging due to the presence of the fermions 
in the fundamental (adjoint) representation. The mass of the eighth pion, related to the 
breaking of the anomaly-free U(l), is given by 

^ = ( 8f + M * P «] ■ (15) 

III. RAISING THE TEMPERATURE 

Next we shall calculate the pressure and the chiral condensates at finite temperature 
up to two loops. The pressure can be calculated following [29 j. At one loop, it equals the 
Bose contribution from a pion gas. Fortunately the two-loop finite temperature results are 
independent on higher derivative terms of the effective zero temperature Lagrangian. To 
this order, therefore, we can be predictive. Beyond two loops one needs the coefficients of 
higher derivative terms that cannot yet be derived from experiments. Lattice simulations, 
in the future, can provide a systematic study of the chiral properties of these theories. 

We summarize below the two-loop results for the pressure for the different patterns of 
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chiral symmetry breaking introduced above: 



Nj - 1 (N 2 - IK 
P(T) = — - — go(m n ,T) gi(m„,T) ; complex (16) 

ft 1\{1AT 1\™2 



(N f + l)(2N f - 1) m ~ 1)(2N; - 1)< 
P(T) = — ^ -go(m nr T) + -* —— 2 gl (m n , T) ; real (17) 

L oIMft n 



(N f - l)(2N f + 1) (Nj ~ 1)(2N/ + \)m\ 



p(T) = go(m n ,T) — — - gi(m„,T) ; pseudoreal(18) 

5 1 

P( T ) = 2^o( m "' r ) + So(m n , T) + -g (m n , T) (19) 

45 wi 2 5 m 2 1 iii 2 

~48 ^ l(m "' 7)2 + 9 7T^( m "' T ) + is -£gifai> T )gi(™n, T) (20) 

+ (8i m ^ + 648 m ^)^^ ; ' (21) 

1 T, 



71 



where 



^= 2 f <22) 



When we have matter in a single representation, the temperature dependent techni- 
quark condensate is found by taking the derivative of the pressure with respect to the 
quark mass. We find 

<QQ> _ 1+ c W 



<QQ>t=o N/FS 9m\ 
where the coefficient c approaches unity as m n goes to zero, and we needed also the zero 

temperature GMOR relation: 

,_ v N f F 2 n m 2 n 

<QQ>r=o = - J -^ JL ■ (24) 
For the first three cases corresponding to a single representation we have: 

(QQ) . ^-Di* Wj-W 

— - — ; = 1 — 7- ; complex (25) 

(QQ) T=0 24N f F 2 n 1152N*F* 

(N f + l)(2N f - 1)T 2 (N 2 f - l)(2N f - 1)T 4 

= 1 - — — — + — r-A ; real (26) 

24N f F 2 n 1152N 2 F\ 

{Nf _ 1)(2N + 1)7 * (N 2 f - l)(2N f + l)T± 

= 1 z — ; pseudoreal (27) 

2AN f F 2 n l\52N 2 Fi v 
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<QQ)(T) 



<QQ>(0) 

1.0 




T (TeV) 



FIG. 1: The temperature dependence of the two-loop chiral condensate for MWT (black solid 
curve), NMWT (red dashed curve), and UMT (blue dotted curve) 



as m n is set to zero. For UMT, we get for the condensate of the fundamental quarks 

<QQ> 



(QQ)t= 



c dp(T) | 4c dp(T) 
+ 2F 2 n dml + 9Fl dml 
2 W2 I 5 



(28) 



m n ,A n =0 I 5 1 F„\ T 



5 Fl 



1 Ft\T 



7 77- (29) 



48 MP Fl \ 1536 2592 f2 5832p4/ f 4 

Let us plot the result for some phenomenologically relevant technicolor models. We 
fix the condensate to the electroweak scale via the relation 



Jbijg F n , with v EW - 246 GeV . 



(30) 



Nf g is the number of flavors gauged under the electroweak symmetry. In Figure [T] we 
plot the condensates as function of the temperature for the MWT, NMWT, and UMT 
models. For (N)MWT we fixed F n to the electroweak scale setting Nf g = 2. For the UMT 
sector only the techniquarks transforming according to the fundamental representation 
of the SU(2) technicolor gauge group are gauged under the electroweak symmetry and 
therefore the decay constant F n of the fundamental sector was fixed to electroweak with 
Nfg = 2. We assumed F n = F n since with two technicolors the fundamental and adjoint 
representations are expected to lead to similar screening effects. 

The temperature effects, to the order we have computed, tend to reduce the size of the 
condensate with respect to the zero temperature value. To estimate the critical temperature 
above which chiral symmetry is expected to restore we set the temperature dependent 
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condensate to zero and read off the associated critical temperature. This procedure is 
not expected to yield a precise estimate of the critical temperature but should capture 
the essential features related to the chiral symmetry breaking pattern on the underlying 
gauge theory 

Therefore, by solving for the temperature where the condensate melts we deduce 
T c = 414, 650, and 480 GeV for MWT, NMWT, and UMT respectively. 

It is interesting also to investigate the general case in which we gauge under the 
electroweak symmetry all doublets: 2N D = Nf = Nf g . The critical temperature can than 
be approximated analytically using the one-loop result once fixed F n to the electroweak 
scale and reads: 

j24N~f V24 

^ - -.F n = — Vew ; complex (31) 



W 1 



J2&N} V24 

-J n = . v EW ; real (32) 



V(N ; + 1)(2N ; - 1) j(N f + l)(2N f -l) 

JlAN~ f V24 

zF n = — Veyj ; pseudoreal (33) 



V(N, - 1)(2N ; + 1) V( N / - W^i + 1) 

It is clear from the formulae just above that the critical temperature is, for large Nf, 
inversely proportional to the number of flavors. 



IV. WALKING CRITICAL TEMPERATURE 



In the previous sections we studied the dependence of the condensate as function of 
the temperature for different technicolor models at fixed number of flavors and matter 
representation. We have also fixed the pion decay constant in order to reproduce the 
correct W gauge boson mass. 

We now deduce the dependence of the critical temperature on the number of flavors 
for a given gauge theory as we approach the critical number of flavors above which the 
theory is expected, at zero temperature, to develop a nonzero infrared fixed point. To 
do so we normalize the walking temperature T c [Nf] to T c [Nf = 2] for the complex and 
pseudoreal representations. We used the two-loop finite temperature result given above 
to deduce the intrinsic dependence on the zero temperature pion decay constant. To get 
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an idea of the Nf dependence we show the one-loop result which fits in one line: 
T c [N f ] fiWf F n [N f ] 



1) 



FJ2] 



T c [2] 

T c [N f ] _ 

T c [2] ^j2(N f - l)(2N f + 1) Fn[2] ' 



complex , 

Fn[Nf] 



pseudoreal . 



(34) 

(35) 
(36) 



From the expressions above it is clear that one needs the explicit dependence of the pion 
decay constant as function of the number flavors. Near the lower end of the conformal 
window one can estimate such a dependence using the Schwinger-Dyson results in the 
rainbow approximation ll33ll35ll4T| . We start with the well known relation: 



2nF n [N f ]* ^Jd[r]J: [N f ] 



(37) 



with d [r] the dimension of the representation of a technifermion species with respect to 
the underlying gauge group and the dynamical mass: 

1 



L [N f ] = A c [N f ]exp 



jaANf] 



1 



(38) 



Since we are comparing a given gauge theory for different number of flavors, in this 
section, we do not indicate the explicit dependence on the fermion representation and 
number of colors which must, however, be taken properly into account when computing 
the walking temperature for each theory 

Here a*[Nf] is the infrared fixed coupling at two loops: 



Pi 

with the following definition of the two-loop beta function: 

p(a) = -fioa 2 - jSia 3 . 



(39) 



We also have: 



7T 



(40) 



(41) 



3C 2 [r] ' 

which is the critical value of the coupling below which one expects the underlying theory 
to develop a nontrivial infrared fixed point and C2I/] is the quadratic casimir. The number 
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a 



0.78 a r 



To 




A, A 



uv 



Log jd 



FIG. 2: The definitions of various scales and couplings in the walking region, with a c close to a*. 



of flavors is assumed to be sufficiently large that the underlying gauge theory breaks 
dynamically chiral symmetry. Since we are approaching the conformal window from 
the broken side «» is the would be infrared fixed point if chiral symmetry would not be 
broken according to the two-loop beta function. A c , following [|4T|, is defined to be the 
energy scale at which the coupling, as function of the renormalization scale, crosses a c . 
Another relevant scale is the one above which the theory is expected to run as in QCD. 
This is the technical scale A uv usually defined such that the coupling evaluated at this 



scale is about 0.78 a c |42|. The above is the standard rainbow approximation scenario and 
the hierarchy of the scales is represented in Fig. |2j 

Because of the existence of a lower number of flavors changing the sign of the second 
coefficient of the two-loop beta function the above approximation of F n as function of 
the number of flavors does not automatically extend to the lowest possible value of 



the number of flavors. When the second coefficient vanishes the exponential in (38) 
approaches unity. Below the associated value of the number flavors we simply take 
£q[N/] = A c [N f ]. 

The explicit dependence of A c on the number of flavors is via the beta function and 
reads: 

mn < ] = a «MLM- <42) 

It is straightforward to see that by approximating the beta function near a*[Nf] leads to: 



A c [N f ] « A 



uv 



a*[Nf]-a c \ft)«<[N/i 



a*[N f ]-0.78a c 



(43) 
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FIG. 3: The normalized critical temperature as a function of Nf for techniquarks in the fundamental 
representation of SU(2) (left) and SU(3) right. 



We can, at this point, plot the normalized walking critical temperature for the two and 
three colors case for fermions in the fundamental representation, respectively, in the left 
and right panels of Fig. |3j 

From the pictures one immediately notices the drop of the temperature near the critical 
number of flavors above which chiral symmetry restores. For very low number of flavors 
we observe the geometric dependence on the number of flavors dictated by chiral pertur- 
bation theory. Our crude approximation in the intermediate region leads to a discontinue 
value of the derivative of the walking critical temperature with the respect to the number 
of flavors at the point when the second coefficient of the beta function vanishes. We 
expect these two distinct regimes, the drop of the critical temperature at low number of 
flavors due to straightforward chiral dynamics, and the one near the lower end of the 
conformal window due to the near conformal dynamics to help disentangling the two 
regions of the phase diagram of gauge theories of fundamental interactions. Within these 
simple approximations for QCD we expect the conformal dynamics to dominate above 
eight flavors while for the two colors case above five flavors [35J. 

The prediction of Braun and Gies J9l on the flavor dependence of the critical tem- 



perature is encoded in our approach in (43), however we observe a further exponential 
suppression related to the dependence of the dynamical mass on the number of flavors in 



(38 1. The main difference is that we have linked the critical temperature to the dynamical 
mass while in the approach of [9] the critical temperature is obtained using dimensional 
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analysis directly from the energy scale in (43). 



V. CONCLUSIONS AND IMPACT ON COSMOLOGY 

We provided the thermal properties of technicolor models using chiral perturbation 
theory At sufficiently large temperatures we explored the melting of the techni-chiral 
condensate. We used the two-loop expressions at nonzero temperature to estimate the 
geometric dependence of the critical temperature on the number of flavors for technimat- 
ter in complex, real and pseudoreal representations of the underlying gauge group and 
also in the presence of multiple distinct matter representations. 

By combining our results with the dependence of the pion decay constant on the 
number of flavors near the lower end of the conformal window, deduced using the 
Schwinger-Dyson approximation, we have shown that the critical temperature walks, i.e. 
decreases near the lower bound of the conformal window. We pointed out the existence 
of two distinct regimes, one at low number of flavors according to which the temperature 
drop as function of flavors is due to chiral pion dynamics and the one for large number 
of flavors which is dominated by the near conformal dynamics. 

We now show that the results obtained in the previous sections are also relevant 
when describing the early universe phase transition via models of dynamical electroweak 
symmetry breaking |6]-l8lH0|. 

In particular the knowledge of the critical temperature, in units of the electroweak scale, 
plays a major role when estimating the order of the phase transition. More precisely one 
estimates the order of the phase transition by computing the ratio, of the thermal value 
of the scalar fermion condensate <fi* (the composite Higgs field) at the bubble nucleation 
temperature T», and T* itself. A strong first order phase transition requires: (j)*/T* > 1. We 
first identify the bubble nucleation temperature with the critical temperature deduced in 
the III section and shown explicitly at one-loop in ( |33] >. 



A first order phase transition is expected to occur © IE! OH for T» < 250 GeV. The 
origin of this result can be easily explained: (p evaluated in the broken phase at the 
minimum decreases with increasing critical temperature and since at zero temperature 
<fi* < <p(T = 0) = 246 GeV the higher is the nucleation (critical) temperature with respect 
to the electroweak scale the smaller is the ratio 0»/T* with respect to unity. 
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Therefore, the phase transition in MWT, NMWT and UMT appears to be either of 
second order or weak first order. 

One can increase the order of the dynamical electroweak phase transition by consid- 



ering the many flavors gauged under the electroweak symmetry In this case from (33) 
it is evident that the critical temperature decreases with the number of flavors favoring 
stronger first order phase transitions. 

Our results provide relevant benchmarks for lattice computations of the conformal 
window ||43Tl84ll at nonzero temperature for any nonsupersymmetric vector-like gauge 
theory with f ermionic matter. There has already been much interest in investigations of the 
thermodynamical properties of the conformal window on the lattice. These studies were 
pioneered by Deuzeman, Lombardo and Pallante [|54ll55| for theories with fermions in the 
fundamental representation and by Kogut and Sinclair [81J for minimal walking models. 
The temperature dependence on the number of flavors stemming from chiral pertrubation 
theory is exact while the model dependence is encoded in the flavor dependence of the 
pion's decay constant. Our predictions can be also used to verify results obtained using 
gauge-gravity duality |85ll86| and, in the future, to further analyze electro-magnetic gauge 
duality E7HS0. 
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